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1.  INTRODUCTION 


Heat  of  explosion  (HEX)  is  usually  defined  as  the  difference  between  the  energies  of  formation  of  a 
material  and  the  enei*’'  :s  of  formation  of  the  explosion  products  (Meyer  1977).  Experimentally,  HEX  is 
measured  as  the  heat  released  when  a  material  is  ignited  and  burned  in  a  bomb  calorimeter  in  an  inert 
atmosphere.  The  measurement  is  made  after  the  products  have  cooled  to  near  room  temperature  by  noting 
t^.e  temperature  rise  of  the  calorimeter  water  jacket;  HEX  is  calculated  from  this  temperature  rise  using 
the  effective  heat  capacity  of  the  calorimeter  body  and  water  jacket.  The  experimentally  determined  HEX 
values  include  contributions  due  to  cooling  of  the  combustion  products  from  the  flame  temperature  (Tf]) 
to  room  temperature  (T^.  and  phase  changes  such  as  water  condensation.  These  contributions  are  not 
generally  available  for  use  in  gun  and  rocket  propulsion,  but,  on  the  other  hand,  HEX  is  much  easier  to 
measure  than  impetus  or  specific  impulse,  and  therefore  is  a  convenient  way  of  ranking  the  energy  content 
of  different  formulations.  Note  that  heat  of  explosion  is  a  misnomer  for  two  reasons:  1)  many  people 
consider  the  term  "heat"  to  imply  constant  pressure  enthalpy  (H);  an  equally  valid  definition  of  this  term 
is  constant  volume  internal  energy  (U).  Since  a  bomb  calorimeter  is  a  constant  volume  device,  the  latter 
usage  of  the  term  will  be  used  in  this  report.  2)  The  sample  usually  undergoes  combustion  rather  than 
ex{dosioa  In  addition,  the  energetic  materials  community  reports  HEX  values  as  positive  numbers,  even 
though  exothermic  energy  release  should  be  negative  according  to  convention  in  the  thermodyr^ics 
cormnunity.  Strictly  speaking  then,  HEX  is  defined  in  this  report  as  the  absolute  value  of  the  total  energy 
released  during  constant  volume  combustion  of  a  material  in  an  inert  atmosphere. 

Measured  HEX  depend  somewhat  on  conditions;  in  particular,  HEX  generally  increases  with  increasing 
loading  density,  or  with  bomb  prepressurization.  To  some  extent,  this  is  predicted  thermodynamically, 
but  there  are  other  potential  contributions:  at  lower  loading  densities  or  prepressurization  levels,  mixing 
with  the  cold  gas  initially  present  in  the  bomb  calorimeter  may  quench  some  of  the  chemistry  before  it 
attains  thermodynamic  equilibrium,  and/or  mixing  with  the  cold  gas  may  raise  the  freeze-out  temperature 
(Tfo)  of  the  cooling  combustion  products.  Tf^  is  the  temperature  during  cooling  at  which  the  ratios  of  the 
various  combustion  products  become  invariant  (Comer  1950),  though  phase  changes  may  still  occur. 

It  is  increasingly  of  interest  to  be  able  to  calculate  HEX  using  a  thennochemical  code  such  as  NASA- 
Lewis  (Gordon  and  McBride  1976;  Svehla  and  McBride  1973)  or  BLAKE  (Freedman  1982).  The  "exact" 
method  for  calculating  HEX  with  a  thermodynamic  code  has  been  described  elsewhere  (Lenchitz  et  al. 
1979)  and  will  be  summarized  below.  Since  the  exact  method  requires  the  Tfo.  which  is  rarely  known. 


1 


and  since,  in  addition,  it  requires  summing  up  the  loom  temperature  energies  of  tlK  products  using  the 
mole  fractions  at  Tf^,  there  is  interest  in  approximate  techniques  which  do  not  require  knowledge  of  Tf^, 
and/or  eliminate  the  need  to  make  use  of  the  composition  of  the  combustion  products.  Three  such 
approximate  techniques  are  described  below.  They  are  illustrated  for  the  case  of  CHNO  chemistry,  but 
may  in  fact  be  especially  useful  for  chemistries  involving  other  elements,  where  essentially  nothing  is 
known  about  Tfo-  In  addition  to  these  three  approximate  techniques,  which  use  general-purpose 
thermochemical  codes,  Baroody  and  Peters  (1990)  have  described  an  approximate  procedure  employing 
a  nockci  motor  code,  where  the  chamber  and  exhaust  conditions  approximate  the  bomb  calorimeter  before 
and  af.er  cooling,  respectively,  and  where  corrections  are  applied  to  compensate  for  the  fact  that  the  rocket 
calculations  are  at  constant  pressure  rather  than  constant  volume. 

HEX  calculations  are  carried  out  on  nitramines,  cyclo-1.3,5-trimethylene-2,4,6-trinitramine  (RDX)  and 
cyclotetramethylene  tetranitramine  (HMX),  as  well  a.  a  mixture  of  nitroglycerine  (NG)  and  nitrocellulose 
(NQ  at  12.6%  nitration  level.  Thermochemical  calculations  were  carried  out  using  CET89,  a  recent 
version  of  the  NASA-Lewis  code  (Gordon  and  McBride  1976;  Svehla  and  McBride  1973)  modified  at  the 
Army  Research  Laboratory  (ARL)  to  run  on  a  personal  computer.  All  calculations  were  run  on  Intel 
80386-  and  80486-based  PCs  operating  under  MS-DOS.  The  procedures  outlined  in  the  following  sections 
would  be  similar  using  other  codes.  In  all  cases,  the  goal  is  to  calculate  or  estimate: 

HEX  .  -  |Up»a  -  u™,| 

where  AUj  is  the  change  in  energy  at  constant  volume  for  the  calorimetric  reaction,  and  Up^  and  U,^ 
are  the  internal  energies  of  the  products  and  reactants,  respectively,  at  room  temperature 
(T,  =  25®C  =  298.15  K).  Analogous  procedures  could  be  used  (using  enthalpies  instead  of  energies)  for 
a  constant  pressure  process.  Table  1  lists  values  for  the  atomic  make-up,  molecular  weight,  and  enthalpy 
of  fonnation  for  the  materials  being  studied. 

2.  ASSIGNMENT  OF  THE  THERMODYNAMIC  QUANTITIES  H  AND  U. 

The  thennodynamic  quantities  H  and  U  are  the  enthalpy  and  internal  energy  state  functions  for  a  given 
species.  These  quantities  are  related  by  the  following  equation: 

H  -  U  +  PV.  (2) 
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Table  1.  Seleaed  Values  for  Materials  Used  in  This  Study 


Material 

C 

H 

N 

0 

MW 

(g/mol) 

AHf° 

(kcal/mol) 

HMX 

D 

8 

8 

8 

296.16 

17.92 

RDX 

6 

6 

6 

222.12 

14.69 

NC,  12.6%  N 

D 

7.550 

2.450 

9.901 

272.40 

-169.20 

NG 

B 

5 

3 

9 

227.09 

-88.60 

where  P  is  pressure  and  V  is  volume.  Equation  2  can  be  modified  to  handle  either  gaseous  or  condensed 
species  as  follows.  For  gaseous  species,  the  ideal  gas  equation  of  state, 

PV  *  nRT,  (3) 

may  be  combined  with  Equation  2  to  give  the  following  equation  relating  H  and  U: 

H  =  U  +  nRT  (gases),  W 

where  R  is  the  gas  constant,  T  is  temperature,  and  n  is  the  number  of  moles  of  gas  involved.  At  room 
temperature,  RT  is  equal  to  592  calAnol.  In  the  case  of  condensed  Giquid  or  solid)  species,  the  PV 
product  of  Equation  2  is  small  enough  that  the  following  approximation  may  be  made: 

H  =  U  (liquids,  solids).  (5) 

Equations  4  and  5  allow  the  assignment  of  internal  energy  in  terms  of  the  corresponding  enthalpy. 

Enthalpies  are  assigned  in  the  codes  by  arbitrarily  assuming  H298.1S  =  where  (A1^98 15 

is  the  standard  heat  of  formation  (relative  to  the  reference  elements  at  298.15  K  and  O.I  MPa  pressure) 
of  the  species  of  interest.  Equation  5  then  becomes 

^  “  ^298.15  ®  ^298.15 

for  liquids  and  solids,  and  Equation  4  becomes  for  gases.  Since  most  energetic  materials  of  interest  exist 

U  =  -  nRT  =  (AHf  )29g.i5  ~  nRT  (7) 
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in  condensed  phases.  Equation  6  is  generaUy  sufficient  to  assign  dK  reactant  energy  (U)  on  the  assigned 
thermodynamic  sctQe.  If  gaseous  product  energies  need  to  be  assigned  (e.g..  see  the  exaa  calculation, 
below),  then  Equation  7  should  be  used. 

A  distinction  needs  to  be  made  between  the  assigned  quantity,  U.  arwl  the  energy  of  fotmadcn,  AUf. 
In  Equations  6  and  7,  U  can  be  thought  of  as  the  internal  energy  of  the  product  of  a  fiormation  reactkm 
(for  the  chemical  species  of  interest)  on  the  assigned  energy  scale;  in  odier  words,  U  represents  an 
"absolute"  energy  on  the  scale  defined  by  setting  H29g,i5  =  (AI^)298.i5'  U  differs  from  die  energy  of 
formation.  AU^  in  that  AUf  is  defined  as  the  difference  of  U  for  die  formadon  species  and  the  U’s  of  the 
elemental  constituents  in  their  standard  states;  Kotlar  (1992)  presents  a  good  discussion  on  the  proper 
interpretation  and  assignment  of  internal  energies.  The  use  of  die  U’s  for  die  reactants  and  products  in 
Equation  1  to  calculate  HEX  is  equivalent  to  the  use  of  the  corresponding  AUf’s  since  the  U  contributions 
(to  AUf)  from  the  elemental  constituents  would  cancel  out  from  the  product  and  reactant  internal  oiergy 
terms.  Often  enthalpies  of  formation  are  more  readily  availatde  than  are  energies  of  formation;  this 
formalism  for  internal  energy  was  chosen  to  simplify  the  calculations  involved  by  usirig  Equatkms  6  and  7. 

3.  HEX  vs.  HEAT  OF  COMBUSTION 

The  oiergy  of  combustion  is  experimentally  measured  through  the  use  of  a  bomb  calorimeter  filled 
with  an  excess  of  oxygen.  For  a  CHNO  species,  die  conbustirm  products  would  be  liquid  H2O  and 
gaseous  CO2  and  N2  according  to  the  fbUowing  equation  (Cox  and  PUdm’,  pp  54): 

C.Hb0^d  y)02(g)  -►  aC02(g)  ♦  .^H2O0)  +  ^N2(g).  (8) 

The  measured  energy  of  omibustion,  AU^,  can  be  converted  to  the  enthalpy  of  combustion,  AH^,  throu^ 
the  use  of  the  following  equation: 

AH^  -  AU^  ♦  AnjRT,  (9) 

where  Arig  is  the  change  in  the  number  of  moles  of  gas  for  the  reaction  and  R  is  the  gas  constant; 
Equation  9  follows  direcdy  from  Equation  4.  When  C,HbO^^  exists  entirely  in  the  solid  or  liquid  state. 
AOg  takes  on  the  value  of  (d/2  +  c/2-  b/4).  The  enthalpy  of  combustion  is  typically  used  to  calculate  the 
endialpy  of  formation  for  the  compound  being  combusted. 
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HEX  is  measured  in  a  similar  manner  as  the  heat  of  combustion  except  that  the  calorimeter  is  filled 
with  an  inert  gas  instead  of  oxygen.  A  simple  calculation  for  estimating  the  enthalpy  of  explosion  (AH^ 
can  be  derived  if  one  assumes  that  the  only  observed  combustion  species  are  CO,  CO2.  H2OO).  and  N2. 
This  assumption  excludes  any  consideration  of  O2.  H2,  or  any  fragmented  species  in  the  combustion 
products;  the  exclusion  of  H2  is  not  unreasonable,  since  the  water-gas  equilibrium  hneeze-out  temperature 
would  typically  be  in  the  1,500-1,800  K  range  (Omellas  1966,  1%8,  1974;  Bemecker  1967),  resulting 
in  equilibrium  constants  which  are  unfavorable  to  the  production  of  H2  (Comer  1950;  Meyer  1977).  The 
explosion  reactitm  is  then  ^proximated  as: 

C^HbO^d  (c  -  a  -  |)C02(g)  +  +  2a  +  |)  CO(g)  +  |  H20(l)  +  ^  N2(g).  (10) 

Derivations  of  Equations  8  and  10  can  be  found  in  Appendix  E. 

The  enthaply  of  combustion  (AH^)  for  Equation  8  can  be  expressed  as: 

AH,  -  aAHf[C02(g)]  +  |  AHf[H20a)]  -  AHflC.HbO^N^]. 

In  a  similar  fashion,  the  enthalpy  of  extdosion  (AH,)  for  Equation  10  can  be  expressed  as: 

AH,  -  (c  -  a  -  i)  AHf[C02(g)]  +  (-c  +  2a  -  .|)AHf[CO(g)] 

+  I  A  Hf  [H20(1)]  -  AHf[C,HbO,Nd]. 

Taking  the  difference  between  Equations  1 1  and  12,  the  following  equation  is  obtained  (solving  for  AH,): 


(11) 


(12) 


AH,  =  AH,  +  (AHf[CO(g)l  -  AHftC02(g)])(2a  +  y  -  c).  (13) 


By  making  the  approximation  HEX  »  -AH,,  Equation  13  can  be  written  (on  a  per  gram  basis)  as: 


HEX 


(-AH,)  -  67.63  •(2a  +  0.5b  -  c) 
MW 


(14) 


where  AH,  is  the  enthalpy  of  combustion  of  the  material  in  kcalAnol,  MW  is  the  molecular  weight  of  the 
compound  in  g/!mol.  and  the  calculated  HEX  is  in  kcal/g.  This  particular  ^rproximation  can  be  justified 
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by  the  fact  that  the  reactions  in  Equations  8  and  10  wouM  be  expected  to  produce  similar  amounts  of 
products  in  the  gas  phase;  the  An^RT  terms  in  Equation  9  would  roughly  cancel  out  if  the  corresponding 
energies  were  used  in  place  of  the  enthalpies  in  Equation  13.  Typically,  it  is  easier  to  find  values  for 
enthalpies  of  formation  than  for  enthalpies  of  combustion.  Where  the  enthalpy  of  combustion  cannot  be 
found  for  a  given  species,  it  may  be  calculated  using  Equation  11.  Table  2  lists  HEX  values  calculated 
using  Equation  14  as  well  as  literature  values  of  HEX  for  HMX,  RDX,  NC,  and  NG.  The  enthalpy  of 
combustion  values  listed  in  Table  2  were  calculated  using  Equation  11  and  are  only  listed  for  use  in 
Equation  14.  It  should  be  noted  that  Equation  14  only  applies  in  the  case  of  pure  compounds  where  the 
heat  of  combustion  can  be  expressed  on  a  molar  basis.  In  the  case  of  a  multi-c(»nponent  formulation. 
Equation  14  could  be  used  to  generate  partial  HEX  for  each  ingredient  which  could  then  be  combined 
on  a  weight  basis  to  estimate  HEX  for  the  formulation  as  a  whole. 


Table  2.  HEX  Values  for  Selected  Energetic  Materials  Calculated  Using  Equation  14 


Material 

-AH, 

(kcalAnol) 

HEX  1 

Eq.  14 
(cal/g) 

Literature 

(cal/g) 

1 

HMX 

667.4 

1,340 

1,356 

Lin  and  Chen  1986 

RDX 

501.8 

1,346 

1,370 

Urbanski  1967 

NC,  12.6%  N 

653.0 

939 

973 

Urbanski  1965 

1 

364.3 

1,753 

1,470 

Urbanski  1965  | 

4.  COMPUTATIONAL  METHODS  USING  THERMOCHEMICAL  CODES 


4. 1  Exact  Calculation  of  HEX.  The  following  equation  is  based  on  the  definition  of  HEX  presented 
in  Equation  I : 


HEX  = 


J-1 


EXjU, 

i-1 _ 

n 


(cal/g). 


(15) 
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where  the  m  product  mole  fractions  (Xj)  in  the  first  summation  are  those  cmnputed  at  the  freeze-out 
temperature  Tfo.  the  room  temperature  internal  eneipes,  Uj,  are  for  the  expected  state  following  cooling 
to  Tf  (e.g.,  H20(1)).  and  division  by  the  product  molecular  weights.  Mj,  gives  HEX  units  of  energy  per 
unit  weight  (e.g.,  cal/g).  For  gaseous  products,  the  Uj  are  related  to  the  enthalpies  of  foimatim  by 
Equation  7;  the  internal  energies  for  use  in  the  first  term  of  Equation  15  are.  therefore,  obtained  by 
subtracting  RT  (592  calAnol)  from  the  room  temperature  enthalpies  of  formation  in  the  JANAF 
thermochemical  tables  (Chase  et  al.  1985).  The  second  term  in  Equation  15  is  simply  the  internal  energy 
(cal/g)  of  the  propellant  or  energetic  material,  obtained  by  summing  the  Uj  for  the  n  ingredients.  For  the 
usual  case  of  a  solid  or  liquid  material.  Equation  6  applies,  so  the  enthalpies  of  formation  can  be  used. 
Most  thermochemical  codes  permit  ingredient  weight  fractions  to  be  input  rather  than  mole  fractions;  most 
codes  output  the  total  energy  of  formation  of  the  formulation,  providing  the  second  term  in  Equation  15 
directly. 

The  exact  procedure  for  calculating  HEX  has  been  described  elsewhere  (Lenchitz  et  al.  1979).  and  is 
briefly  summarized  as  follows:  If  the  code  supports  an  assigned  temperature  and  volume  mode,  as  does 
the  NASA-Lewis  code,  specify  a  volume  (or  loading  density,  p)  and  T,  where  an  assumed  Tf^  is  used  for 
T.  and  run  the  code  to  produce  the  minimum  energy  (equilibrium)  combustion  product  composition.  Fbr 
a  CHNO  formulation,  the  only  significant  products  are  usually  CO.  CO2,  H2,  H20(g4),  N2,  CH4,  NH3, 
and  C(gr);  where  the  methane,  ammonia,  and  grr^te  result  from  oxygen  deficiency  in  the  energetic 
material  or  formulation.  The  predicted  mole  fractions  at  Tf^  are  then  used,  together  with  the  molecular 
weights  and  room  temperature  internal  energies  to  perform  the  summation  in  the  first  term  of  Equation  15. 
In  the  original  implementation  (Lenchitz  et  al.  1979),  this  was  done  by  taking  the  produa  mole  fractions 
at  Tfo  from  the  thermochemical  code  output,  and  inputting  them  into  a  small  program  running  on  a 
programmable  calculator.  More  recently,  we  have  done  the  summation  by  downloading  the 
thermochemical  code  output  into  a  spreadsheet  (see  Appendix  A).  Yet  another  procedure,  for  some 
thermochemical  codes,  is  to  enter  the  gas  composition  returned  from  the  assigned  temperature  and  volume 
calculation  as  input  to  an  assigned  energy  and  volume  calculation,  which  is  then  run  at  T^  =  298.15  K; 
when  specifying  the  composition,  any  gaseous  H2O  should  be  condensed  to  liquid  for  the  room 
temperature  calculation.  The  final  composition  returned  from  the  assigned  energy  and  volume  calculation 
is  ignored,  only  making  use  of  the  output  U  for  Uk  mixture;  since  the  calculation  is  run  specifying 
constant  U,  the  calculation  is  essentially  a  lazy  way  of  summing  up  the  internal  energy  of  the  ir^t  (Tj^ 
composition  at  room  temperature.  Of  course,  the  ideal  procedure  is  to  modify  the  thermodiemical  code 
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to  calculate  a  HEX  automatically  for  oi^  or  more  ^cified  fireeze-out  temperatures.  This  modification 
has  been  implemented  in  a  recent  version  of  the  BLAKE  thermochemical  code. 


Table  3  shows  calculated  HEX  for  HMX  and  for  a  60/40  mixture  (by  weight)  of  nitrocellulose  (12.6% 
N)AuttDglycerine  (NC/NG)  using  the  exact  procedure.  The  variation  of  HEX  with  loading  density,  p,  is 
readily  tyqrarent  in  the  data  of  Table  3,  and  is  due  to  increased  formation  of  C(gr),  CH4,  and  NH3  with 
decreasing  freeze-out  temperature. 


If  a  code  does  not  support  an  assigned  temperature  and  volume  (or  p)  mode,  the  problem  then 
becomes  one  of  estimating  the  pressure  at  Tf^  (i.e.,  the  freeze-out  pressure,  Pf^)  so  that  an  assigned 
temperature  and  pressure  calculation  can  then  be  carried  out  to  find  the  composition  at  freeze  out  The 
Pfo  can  be  estimated  from  (Lenchitz  et  al.  1979): 


(16) 


where  Tf]  and  Pf]  are  the  adiabatic  flame  temperature  and  pressure  (i.e.,  the  T  and  P  before  any  cooling 
occurs).  Tfi  and  Pfi  are  computed  using  an  assigned  energy  and  volume  calculation.  Equation  16 
overestimates  Pfo  somewhat,  since,  in  general,  condensed  species  form  at  Tf^,  and/or  the  average  molecular 
weight  of  the  gases  increases  on  cooling  toTfo.  This  leads  to  the  resulting  HEX  being  a  little  too  high; 
for  typical  CHNO  propellants  the  error  is  usually  only  20  cal/g  or  less  (Lenchitz  et  al.  1979).  This  error 
can  be  eliminated  by  using  an  iterative  procedure  in  place  of  Equation  16.  This  is  done  as  follows.  Hrst, 
estimate  Pfo  from  the  ideal  gas  equation  of  state  P  =  pRT/M  (for  P  in  atm  and  p  in  g/cm^, 
R  =  82.06  cm^-atm/hiol-K)  for  the  p  and  T^^  of  interest,  using  a  reasonable  guess  for  the  molecular  weight 
M.  Run  the  thermodynamic  code  at  assigned  temperature  and  pressure  using  this  P  and  T.  The  resulting 
returned  M  is  then  used  with  the  ideal  gas  equation  of  state  to  obtain  an  improved  estimate  of  Pf^,  this 
procedure  being  repeated  until  the  computed  p  is  within,  say,  1%  of  the  desired  p.  If  the  initial  estimate 
of  M  is  close,  generally  only  two  runs  will  be  required;  otherwise  three  should  be  sufficient  to  produce 
HEX  identical  to  those  with  a  code  having  the  assigned  temperature  and  volume  option. 


4.2  The  CydT  Aimroximate  Method.  The  first  sq^roximate  method  we  will  describe  preserves  the 
concept  of  a  freeze-out  temperature,  while  eliminating  the  need  to  sum  up  rotnn  temperature  internal 
energies  for  the  product  species  in  equilibrium  at  Starting  hom  Equation  1,  we  have: 
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HEX 


|(Ufc  -  c,at)  -  u^|. 


(17) 


-  u 


react 


where  and  Uf^  are  the  specific  heat  and  total  internal  eneigy  of  die  product  gases  at  Tfo.  respectively. 


Table  3.  HEX  (cal/g).  Exact  Calculation 


HMX* 

Tfo(K) 

NC/NG** 

Tfo(K) 

p  (g/cm^ 

900 

1,100 

1,300 

1,500 

900 

1,100 

1,300 

1,500 

0.01 

1,504 

1,382 

1,356 

1,356 

1,408 

1,291 

1,278 

1,278 

1,540 

1,445 

1,372 

1,357 

1,447 

1,341 

1,284 

1,279 

1,549 

1,469 

1,393 

1,361 

1,457 

1,365 

1,295 

1,280 

1  0.20 

1,556 

1,491 

1,419 

1,371 

1,464 

1,391 

1,315 

1,283 

1  0.50 

1,562 

1,517 

1,454 

1,398 

1,470 

1,421 

1,347 

1,297 

■  HMX:  U  -  AH/*  >  +  17.92  kcaVmol  -  4A0.S2  cal/g. 

•*  NC  (1Z6«.N);  U  •  AH/* »  -169.2  kcal/mol  -  -621.0  cal/g. 

NG:  U  •  AH/*  >  -88.60  kcal/hiol  >  390.2  cal/g. 

For  60%  NC  40%  NG.  U  «  -528.7  cal/g. 

and  AT  =  Tfo  -  298.15.  Here,  we  are  essentially  approximating  the  internal  energy  of  the  products  after 
cooling  to  room  temperature  by  subtracting  C^AT  hnm  the  internal  energy  output  by  the  thermochemical 
code  for  the  products  at  the  heeze-out  temperature:  U{^  in  Equation  17  is  calculated  firom  an  assigned 
temperature  and  volume  calculation  run  at  a  specified  p  and  by  setting  T  =  Tfo.  Since  specific  heat  is 
actually  a  temperature-dependent  quantity,  an  exact  procedure  would  involve  evaluation  of  JCydT  radier 
than  CyAT.  Also,  this  approximation  does  not  include  the  energy  of  phase  changes  such  as  water 
condensatioa  For  typical  CHNO  formulations,  the  net  result  of  these  tqiproximations  is  that  die  HEX  are 
about  40  cal/g  lower  than  for  the  exaa  method,  as  shown  in  TaUe  4.  Fbr  HMX  and  NC/NG  over  a  range 
of  loading  densities,  simply  adding  40  cal/g  to  the  HEX  from  Equation  17  gives  values  that  ate  within  a 
fraction  of  a  percent  of  those  from  the  exact  method  described  above. 

Depending  on  the  dietmochemical  code,  some  caution  is  required  in  determining  the  iq)|Ht^ate  Cy 
to  use  in  Equation  17.  Many  codes  output  only  C^,  even  for  constant  volume  calculations.  In  this  case, 
Cy  can  be  estimated  from: 
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Table  4.  CydT  Method  Approximate  HEX  (cal/g)  and  Percentage  Difference  from  "Exact"  Calculatkm 
(Table  3) 


Tfo 

(K) 

p  (g/cm^ 

900 

1,100 

1,300 

1,500 

HMX 

0.01 

1,452  (-3.42) 

1,341  (-2.94) 

1,317  (-2.90) 

1,321  (-2J3) 

0.05 

1,483  (-3.72) 

1,406  (-2.64) 

1,333  (-2.82) 

1,323  (-2J3) 

0.10 

1,491  (-3.79) 

1,433  (-2.49) 

1,356  (-2.66) 

1.327  (-2.48) 

0.20 

1,496  (-3.84) 

1.455  (-2.39) 

1,384  (-2.44) 

1,338  (-2.39) 

0.50 

1,501  (-3.89) 

1,479  (-2.52) 

1,423  (-2.12) 

1,368  (-2.12) 

1  HMX  -1-  40  cal/g 

0.01 

1,492  (-0.76) 

1,381  (-0.05) 

1,357  (0.04) 

1,361  (0.42) 

0.05 

1,523  (-1.12) 

1,446  (0.13) 

1.373  (0.04) 

1,363  (0.42) 

0.10 

1.531  (-1.21) 

1,473  (0.23) 

1,396  (0.22) 

1.367  (0.46) 

0.20 

1,536  (-1.27) 

1,495  (0.29) 

1,424  (0.37) 

1,378  (0.52) 

0.50 

1,541  (-1.33) 

1,519  (0.12) 

1,463  (0.63) 

1,408  (0.74) 

1  NC/NG 

0.01 

1,353  (-3.89) 

1,243  (-3.70) 

1231  (-3.68) 

1237  (-324) 

0.05 

1,386  (-421) 

1,297  (-322) 

1238  (-3.60) 

1237  (-323) 

0.10 

1,394  (-4.30) 

1,325  (-2.95) 

1251  (-3.45) 

1238  (-323) 

020 

1,400  (-4.36) 

1,351  (-2.89) 

1273  (-3.17) 

1243  (-3.16) 

0.50 

1.405  (-4.42) 

1,378  (-3.02) 

1,310  (-2.70) 

1259  (-2.92) 

NC/NG,  +  40  cal/g 

0.01 

1,393  (-1.05) 

1,243  (-0.60) 

1271  (-0.55) 

1,277  (-0.11) 

0.05 

1,426  (-1.45) 

1,337  (-023) 

1,278  (-0.49) 

1277  (-0.10) 

0.10 

1,434  (-1.55) 

1,365  (-0.02) 

1,291  (-0.36) 

1278  (-0.10) 

0.20 

1,440  (-1.63) 

1,391  (-0.01) 

1,313  (-0.13) 

1283  (-0.04) 

0.50 

1,445  (-1.70) 

1,418  (-0.20) 

1,350  (0.27) 

1299  (-0.17) 
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(18) 


Cv 


M 


where  M  is  the  average  molecular  weight  of  the  products.  In  the  case  of  the  NASA-Lewis  code  ou^t, 
Cy  can  also  be  determined  from: 


Cv 


(19) 


(dlnV/dlnP>].  is  the  returned  compressibility  and  is  die  isentropic  gamma,  which  is  goierally  within  a 
few  percent  of  y  =  Cp/C^.  With  the  NASA-Lewis  code,  tlrere  are,  however,  additional  complications: 
the  Cp  and  y^  returned  by  the  code  are  equilibrium  rather  than  frozen  values,  in  that  they  include 
contributions  due  to  shifting  chemical  equilibria  in  the  combustion  products.  For  propellants  and 
exfdosives,  it  has  been  shown  (Fifer  and  Cohen  1977)  that  for  sufficioitly  hi^  pressures  and  low  loading 
densities,  the  equilibrium  Cp  may  be  consideraUy  higher  than  the  frozen  C^  needed  for  the  inesent 
calculation.  The  use  of  equilibrium  rather  than  ftozoi  properties  will  result  in  anomalously  hi^  HEX. 
The  safest  way  to  preclude  this  potential  problem  iKlien  using  NASA-Lewis  is  to  request  that  tran^rt 
values  also  be  returned;  this  results  in  frozm  specific  heats  being  returned,  resulting  in  HEX  values  similar 
to  those  determined  using  codes  that  output  only  ftozen  parameters.  In  this  case,  only  Equation  18  should 
be  used  to  calculate  a  frozen  Cp;  under  frozen  conditions.  Equation  19  reduces  to  C^  =  Cp/y,  Mriiere  y  is 
not  explicitly  calculated  by  the  code.  Note,  however,  that  for  systems  which  produce  condensed  products, 
the  returned  transport  properties  (including  C^)  are  for  the  gas  phase  only;  in  addition,  the  M  returned  by 
the  NASA-Lewis  code  (for  use  in  Equation  18)  is  not  a  true  average  molecular  weight,  but  rather 
corresponds  to  a  quantity  defined  by  total  weight  (both  gaseous  and  condensed)  divided  by  the  moles  of 
gas  only.  A  sample  calculation  for  this  method  whidi  makes  use  of  Equation  18  is  presented  in 
Appendix  B  of  this  report.  When  dealing  with  systems  that  produce  significant  amounts  of  conda)s«l 
products,  it  is  recommended  to  either  use  a  code  (other  than  NASA-Lewis)  which  returns  true  average 
molecular  weights  or  use  a  different  approximate  method. 


4.3  The  Room  Temperature  Method.  This  is  the  first  of  two  methods  that  not  only  do  not  require 
weighted  summation  of  product  internal  energies  (Equation  IS),  but,  in  addition,  do  not  require  freeze-out 
temperatures  to  be  specified.  Such  approximate  methods  are  particularly  approfniate  for  diemistries  or 
calorimetric  conditions  where  there  is  little  informatimi  about  the  freeze-out  temperature(s). 
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In  the  room  temperature  {^proximate  method,  it  is  assumed  that  there  is  ix>  fteeze-out  of  chemistry 
on  cooling  of  the  produc  '  to  This  is  equivalent  to  saying  that  =  298.15  K.  To  do  the 

calculation,  an  assigned  temperature  and  volume  (or  density)  calculation  is  carried  out  for  the  desired 
loading  densities,  specifying  T  =  298.15.  From  the  returned  for  the  combustion  products  at  298.15 
K,  is  subtracted  (Equation  1)  to  give  HEX  (see  Appendix  C  for  an  example  calculation).  In  using 
this  method,  we  usually  restricted  the  products  to  CO.  CO2.  H2.  N2.  H2O,  and  H20(l),  since  otlrerwise  at 
this  low  temperature  unrealistically  large  amounts  of  C(gr),  CH4.  NH3,  and  other  products  are  predicted. 
As  shown  in  Table  5,  it  is  then  found  that  the  computed  HEX  are  virtually  independent  of  loading  density 
for  CHNO  systems,  with  HEX  =  1,355  cal/g  for  HMX  and  1,278  cal/g  for  NC/NG.  The  HEX  are  similar 
to  those  calculated  in  the  exact  method  with  a  high  (e.g.,  1400  K)  rather  than  low  (e.g.,  9(X)  K)  Tfo-  This 
results  primarily  from  the  exclusion  of  C(gr).  CH4,  and  NH3.  and  to  a  much  smaller  extent,  from  alteration 
of  the  water-gas  equilibrium  (CO  +  H2O  ••  CO2  +  H2)  due  to  fonnation  of  H20(l)  at  equilibrium. 
Allowing  all  species,  the  HEX  increases  to  about  1 .630  cal/g  for  HMX,  and  about  1 ,524  cal/g  for  NC/NG. 
As  in  the  exact  and  CydT  approximate  methods,  a  dummy  assigned  energy  and  volume  calculation  can 
be  used,  if  desired,  to  calculate  from  the  energies  of  formation  for  the  separate  ingredients  and  an 
assigned  temperature  and  pressure  calculation  can  be  used,  with  estimated  P.  to  obtain  the  product  Up^ 
if  the  code  does  not  offer  the  assigned  temperature  and  volume  mode.  For  CHNO  chemistry  with  CH4, 
NH3,  and  C(gr)  excluded,  the  density-  (pressure-)  independence  of  the  computed  Up^  means  that  an 
accurate  estimate  of  P  is  not  necessary,  so  no  iterations  are  required. 

4.4  The  Inert  Diluent  Method.  The  other  q^roximate  method  that  avoids  both  the  summation  of 
Equation  15  and  the  need  to  specify  Tf^  is  the  inert  diluent  method.  In  this  method,  a  large  excess  of  inert 
gas  (e.g..  Argon)  is  added  to  the  energetic  material,  and  an  assigned  energy  and  volume  (or  density) 
calculation  run.  Tire  effect  of  the  diluent  is  to  provide  a  heat  sink  for  the  reaction  exothermicity,  resulting 
in  computed  adiabatic  flame  temperatures  that  are  only  a  few  degrees  above  room  temperature.  HEX  is 
determined  using  the  following  equation; 

HEX  -  .  N/c.dT  .  (20) 

where  AT  is  the  computed  temperature  rise  (Tf,  -  Tj„j,j,|),  N  is  the  reciprocal  of  the  weight  fraction  of 
propellant  or  explosive  in  the  mixture,  and  Cp  and  y  are  as  have  been  previously  defirred.  The  smaU  AT 
serves  three  purposes:  a)  it  allows  condensables  to  condense,  as  in  the  calorimeter,  b)  it  reduces  the  error 
associated  with  assuming  a  constant  specific  heat;  and  c)  for  the  NASA-Lewis  code,  the  low  T^]  makes 
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it  less  likely  that  die  Cp  will  be  too  large  due  to  equilibrium  reaction  effects,  although  in  some  cases, 
phase  equilibria  could  result  in  equilibrium  being  larger  than  the  frozen  ones.  Unlike  all  of  the  other 
methods  described  above,  this  method  does  not  require  the  explicit  subtraction  of  from  an 
afqxoximation  to  Up^,  but  instead  evaluates  HEX  according  to  the  right-hand  side  of  Equation  20.  An 
examfde  calculation  for  this  method  is  presented  in  Ai^ndix  D.  As  for  the  room  temperature  method, 
it  is  necessary  to  exclude  many  species  that  would  otherwise  be  predicted  in  unreasonable  amounts, 
including  CH4,  NH3,  and  C(gr),  and,  as  for  the  room  temperature  method,  this  results  in  HEX  values  that 
are  almost  totally  independent  of  loading  density. 


Table  S.  Summary  of  All  Methods  for  HEX  (cal/g)  at  Selected  Conditions 


Tfo(K) 

"Exact” 

"CvdT  +  40" 

"Room  T' 

No.  Species 

"Inert  Diluent" 

No.  Species 

1  p(g/cm^) 

1400  K 

900K 

1400  K 

6 

all 

6  +  At 

aU 

Equation 

I4W 

1  hmx  I 

0.01 

1A04 

1456 

1,492 

1461 

1454 

1.629 

1455 

1.495 

1440 

1456 

0.0s 

1A40 

1457 

1423 

1463 

1455 

1,630 

1455 

1.495 

1440 

1456 

0.10 

1A49 

1461 

1431 

1467 

1455 

1,630 

1455 

1.495 

1440 

1456 

0.20 

1455 

1471 

1436 

1478 

1455 

1,630 

1455 

1.495 

1440 

1456 

OAO 

1462 

1498 

1441 

1455 

1,630 

1455 

1.632* 

1440 

1456 

1  NC/NG  1 

0.01 

BS9 

U78 

1493 

1477 

1476 

1422 

1471 

1497 

1465 

— 

0.0s 

1.447 

1479 

1,426 

1477 

1477 

1423 

1471 

1497 

1465 

— 

0.10 

1.457 

1480 

1,434 

1478 

1478 

1424 

1471 

1497 

1465 

0.20 

1,464 

1483 

1,440 

1483 

1478 

1424 

1471 

1497 

1465 

— 

0.50 

1,470 

1497 

1,445 

1499 

1478 

1424 

1471 

1419* 

1465 

— 

*  H2O  0)  appean  al  this  and  higher  p. 

^  From  Table  2.  Values  for  N(^G  are  calculated  by  taking  60%  HEXur  40%  HEXjijO  calculated  using  Equation  14. 

HEXh(.=  I,7S3  cal/g. 

*  No  loadiiig  density  wu  given  for  the  literature  value  for  HMX  (Table  2). 

^  It  is  not  intended  to  imply  that  these  values  are  independent  of  loading  density. 

A  convenient  diluent  is  Argon.  Using  1  part  formulation  +  999  parts  (by  weight)  of  Argon 
(N  =  KXX)),  the  AT  is  typically  20-22  K.  For  sufficient  accuracy,  it  is  necessary  to  have  the  code  ouqxit 
Tfi  to  ±0.1,  or  preferably  ±0.01  K.  Less  diluent  (e.g.,  N  =  100)  can  be  used  if  the  source  code  is  not 
available  to  make  this  format  modification.  Using  999  parts  of  Argon,  and  allowing  only  CO,  CO2.  H2. 
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N2,  HjO,  and  HjOO)  to  form.  Equation  20  gives  HEX  values  of  1,355  cal/g  for  HMX,  and  1,271  cal/g 
for  NCVNG.  For  HMX,  as  shown  in  Table  5,  this  value  is  identical  to  that  using  the  room  temperature 
method;  for  NC/NG,  it  is  seven  cal/g  less.  If  all  species  are  allowed,  HEX  values  for  HMX  are  about 
140  cal/g  higher,  and  for  NC/NG  about  126  cal/g  higher,  but  unrealistically  high  amounts  of  graj^te 
(e.g.,  20-30  mol  %)  are  predicted.  With  all  species  allowed  and  999  parts  Argon,  there  is  a  sudden 
additional  increase  of  122  O^C/NG)  to  137  (HMX)  cal/g  in  going  from  p  =  0.2  to  p  =  0.5  g/cm^;  this  is 
due  to  the  appearance  of  H2OG)  at  this  loading  density.  Of  course,  H20(l)  would  appear  at  much  lower 
loading  densities  if  it  were  not  for  the  Argon  diluent.  This  points  out  the  value  of  running  the  code  for 
extremely  high  loading  densities  (e.g.,  p  =  10  or  100  g/cm^);  although  these  are  unachievably  high  loading 
densities,  they  force  the  water  to  condense  as  it  would  in  the  absence  of  the  inert  diluent. 

With  the  inert  diluent  method,  it  is  also  possible  to  shift  the  ratios  of  predicted  combustion  products 
by  using  one  of  the  combustion  products  as  the  (not  exactly  inert)  diluent.  For  example,  for  HMX,  using 
CO2  as  the  diluent  increases  CO  somewhat,  lowering  HEX  7  cal/g  compared  to  Argon,  using  CO  as 
diluent  (N  =  700)  suppresses  all  water  formation  (HEX  unchanged),  using  H2O  as  diluent  (N  =  500, 
H2OO)  not  allowed)  suppresses  all  CO  (HEX  increased  1  cal/g),  and  using  H2  as  diluent  (N  =  50  due  to 
much  lower  Cp  for  H2)  causes  more  H20(l)  to  form  and  increases  CO/CO2,  causing  HEX  to  decrease  by 
12  to  28  cal/g  depending  on  p.  While  these  changes  are  small  due  to  the  low  AU  for  the  water  gas 
(CO  +  H2O  «■  CO2  +  H2)  reaction,  for  other  chemistries,  use  of  a  combustion  product  as  diluent  may  in 
some  cases  help  force  a  better  match  between  predicted  products  and  those  recovered  from  the  calorimeter, 
producing  a  better  match  between  experimental  and  theoretical  HEX. 

5.  SUMMARY 

Three  approximate  methods  to  estimate  HEX  for  energetic  materials  using  a  standard  thermochemical 
code  have  been  outlined.  The  results  of  these  methods  have  been  compared  to  those  obtained  using  the 
exaa  method  for  calculating  HEX  at  a  given  freeze-out  temperature.  One  of  these  methods  requires  the 
use  of  a  freeze-out  temperature  while  the  other  two  do  not  The  three  approximate  methods  and  the  exact 
method  are  reviewed  below. 

The  exact  method  is  based  on  running  an  assigned  temperature  and  volume  calculation  at  an  assumed 
freeze-out  temperature.  The  calculated  product  composition  (at  T^  is  used  to  calculate  a  rotMn 
temperature  internal  energy  based  on  the  room  temperature  internal  energies  for  the  product  species.  Any 
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gaseous  water  is  converted  to  liquid  water  prior  to  calculating  the  internal  energy  of  the  combustion 
products.  HEX  is  calculated  by  subtracting  the  internal  energy  of  the  uncombusted  material  from  that  of 
its  combustion  products. 

The  CydT  approximate  method  is  also  based  on  running  an  assigned  temperature  and  volume 
calculation  at  an  assumed  freeze-out  temperature.  This  method  avoids  the  summation  of  combustion 
product  internal  energies  by  subtracting  the  quantity  C^AT  from  the  calculated  internal  energy  at  the 
freeze-out  temperature  to  get  the  approximate  internal  energy  of  the  combustion  products  at  room 
temperature.  As  in  the  exact  method.  HEX  is  calculated  by  subtracting  the  internal  energy  of  the 
uncombusted  material  from  that  of  its  combustion  products.  The  AT  is  Tf^  -  298.15  K  and  the  is  the 
constant  volume  heat  capacity,  either  returned  by  the  code  or  calculated  as  a  derived  quantity  from  other 
values  returned  by  the  code.  One  needs  to  be  careful  to  choose  a  code  that  can  handle  condensed  species 
if  they  are  produced.  For  the  CHNO  chemistries  looked  at  in  this  report,  HEX  values  calculated  were 
about  40  cal/g  lower  than  those  using  the  exact  method. 

The  room  temperature  method  requires  no  preassigned  freeze-out  temperature,  and  like  the  CydT 
ai^rroximate  method,  avoids  the  summation  of  comlmstion  product  internal  energies.  As  is  the  case  with 
the  exact  method  and  the  C^dT  approximate  method,  this  method  is  based  on  an  assigned  temperamre  and 
volume  calculation,  but  unlike  these  other  two  methods,  the  calculation  is  run  at  298.15  K.  The 
combustion  produa  internal  energy  is  calculated  directly  from  the  code,  and  HEX  is  determined  by 
subtracting  the  internal  energy  of  the  uncombusted  material  from  that  of  its  combustion  products.  There 
is  a  need  to  restrict  the  permissable  product  species  when  using  this  method  to  prevent  excess  formadtHi 
of  C(gr),  CH^,  NH3.  and  similar  combustion  products  in  the  code.  The  calculated  HEX  is  virtually 
independent  of  loading  density  when  this  method  is  used. 

The  inert  diluent  method,  like  the  room  temperature  method  requires  no  assumed  fteeze-out 
temperature  and  also  avoids  the  summation  of  combustion  product  internal  energies.  In  this  method,  a 
large  excess  of  inert  gas  (e.g..  Argon)  is  added  to  the  energetic  material  and  an  assigned  energy  and 
volume  calculation  is  run.  The  eflea  of  the  diluent  is  to  limit  the  temperature  rise  to  a  few  tens  of 
degrees.  HEX  is  determined  using  a  few  quantities  repotted  by  the  code  without  the  need  to  subtract 
^react  ^  approximation  to  Up^.  As  with  the  room  temperature  method,  there  is  a  need  to  omit 
species  such  as  C(gr),  CH4,  and  NH3,  and  the  calculated  HEX  values  are  almost  totally  indepeiKlent  of 
loading  doisity.  Due  to  the  small  temperature  rise,  the  code  is  required  to  output  temperature  to  ±0.1  K 
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or  ±0.01  K.  Water  can  be  forced  to  condense  out  using  large  loading  densities  of  10  or  100  g/cn?,  but 
since  HEX  is  relatively  unaffected  by  loading  density,  this  does  not  make  a  large  difference. 

Since  two  of  the  approximate  techniques  (room  temperature  and  inert  diluent)  eliminate  the  need  to 
specify  a  composition  freeze-out  temperature,  these  may  be  appropriate  for  situations  where  ^fo 
unknown,  especially  in  cases  which  cannot  be  described  by  simple  CHNO  chemistry  (such  as  metallized 
propellants  or  pyrotechnics).  It  is  interesting  to  note  that  these  two  methods  give  similar  values  which 
seem  to  be  quite  independent  of  loading  density  if  the  number  of  produa  species  allowed  is  kept  small. 
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APPENDIX  A: 

EXAMPLE  OF  THE  EXACT  METHOD  FOR  HMX,  T^,  *  900  K. 
LOADING  DENSITY  OF  0.1  %/ca? 
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INTENTIONALLY  LEFT  BLANK. 
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Stan  by  miming  an  assigned  temperature  and  volume  (TV  =  T)  calculation  at  the  desired  freeze-out 
temperature  and  loading  density.  The  input  file  looks  as  follows: 

REACTANTS 

C  4.  H  8.  N  8.  08.  .  1.0  17924. 5S  298.15  HMX 

NAMELISTS 

&INPT2  TV=T, RHO-0 . 1, T-900, TRACE-1 .E-03/ 

The  following  calculated  values  are  extracted  from  the  output  file: 


THERMODYNAMIC  PROPERTIES 


P,  ATM 

T,  DEG  K 
RHO,  G/CC 
H,  CAL/G 

U,  CAL/G 
G,  CAL/G 

S,  CAL/{G){K) 

M,  MOL  WT 
(DLV/DLP)T 
(DLV/DLT)P 
CP,  CAL/  (G)  (K) 
GAMMA  (S) 

SON  VEL, M/SEC 
IMPETUS  (J/G) 

MOLE  FRACTIONS 

CH4 

CO 

C02 

H2 

H20 

NH3 

N2 

C(GR) 


254.55 

900.00 

1.0000-1 

-1198.13 

-1259.77 

-2695.48 

1.6637 

29.013 

-1.01010 

1.1610 

.5522 

1.1864 

553.2 

257.9 


6.0821-2  <“ 
1.2956-2  <— 

2.5527-1  <~ 

2.6358-2  <~ 

2.2487-1  <— 

8.5717-4  <— 

3.7375-1  <“ 

4.5115-2  <— 


The  calculated  mole  fractions,  heats  of  formation,  and  molecular  weights  are  loaded  into  a  si»eadsheet  to 
calculate  the  internal  energy  (at  298.15  K)  of  the  comlmstion  products  for  the  equilibrium  composition 
calculated  at  the  freeze-out  temperature.  All  water  is  forced  to  condoise  out  (the  heat  of  foimation  for 
liquid  water  should  be  used). 
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Species 

Heat  of 
Formation 
(kJ/mole) 

Heat  of 
Formation 
(kcal/mole) 

-RT 

U<298.15  K) 
(kcal/mole) 

CH4 

-74.873 

-17.895 

-0.592 

-18.487 

CO 

-110.527 

-26.417 

-0.592 

-27.009 

C02 

-393.522 

-94.054 

-0.592 

-94.646 

H2 

0.000 

0.000 

-0.592 

-0.592 

H20(L) 

-285.830 

-68.315 

-68.315 

NH3 

-45.898 

-10.970 

-0.592 

-11.562 

N2 

0.000 

0.000 

-0.592 

-0.592 

C<GR) 

0.000 

0.000 

0.000 

Species 

X 

M 

U 

(cal/mole) 

X*U 

X*M 

CH4 

6.0821E-02 

16.043 

-18487 

-1124.4 

0.9757 

CO 

1.2956E-02 

28.010 

-27009 

-349.9 

0.3629 

C02 

2.5527E-01 

44.010 

-94646 

-24160.3 

11.2344 

H2 

2.6358E-02 

2.016 

-592 

-15.6 

0.0531 

H20(L) 

2.2487E-01 

18.015 

-68315 

-15362.0 

4.0511 

NH3 

8.5717E-04 

17.031 

-11562 

-9.9 

0.0146 

N2 

3.7375E-01 

28.013 

-592 

-221.3 

10.4700 

C(GR) 

4.5115E-02 

12.011 

0 

0.0 

0.5419 

Summations 

-41243.4 

27.7037 

Sum(X*U) /Suin(X*M)  -1488.73  cal/g 


The  calculated  value  of  -1,488.73  cal/g  is  the  internal  energy  of  the  combustion  products  at  298.1S  K 
(using  the  fiteze-out  composition).  The  internal  energy  of  the  uncombusted  HMX  is  subtracted  from  this 
value  to  determine  HEX: 

HEX  -  1-1488.73  -  +  60.521  -  1549  cal/g 
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APPENDIX  B: 

EXAMPLE  OF  C^dT  APPROXIMATE  MEITIOD  FOR  HMX,  *  900  K, 
LOADING  DENSITY  OF  0.1  g/cm^ 
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IKTENTIONALLY  LEFT  BLANK. 
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Stait  by  running  an  assigned  temperature  and  volume  (TV  =  T)  calculation  with  transport  properties 
(TRNSFT  =  T)  at  the  desired  freeze-out  temperature  (900  K)  and  loading  density  (0.1  g/on^).  The  ii^ 
file  looks  as  follows: 


REACTANTS 

C  4.  H  8.  N  8.  08.  .  1.0  17924. 5S  298.15  HMX 

NAMELISTS 

& INPT2  TV-T , TRNSPT-T , T- 9  0  0 , RHO-0 . 1 , TRACE-1 . E-0  3 / 

The  following  calculated  values  are  extracted  from  the  output  file: 


THERMODYNAMIC 

PROPERTIES 

P,  ATM 

254.55 

T,  DEG  K 

900.00 

RHO,  G/CC 

1.0000-1 

H,  CAL/G 

-1198.13 

U,  CAL/G 

-1259.77 

G,  CAL/G 

-2695.48 

S,  CAL/ (G) (K) 

1.6637 

M,  MOL  WT 

29.013 

(DLV/DLP)T 

-1.01010 

(DLV/DLT) P 

1.1610 

CP,  CAL/(G)(K) 

.5522 

GAMMA  (S) 

1.1864 

SON  VEL,M/SEC 

553,2 

IMPETUS  (J/G) 

257.9 

TRANSPORT  PROPERTIES  (GASES  ONLY) 

CONDUCTIVITY  IN  UNITS  OF  MILLICALORIES/ (CM) (K) (SEC) 

VISC,MILLIPOISE  .37253 

WITH  EQUILIBRIUM  REACTIONS 

CP,  CAL/(G)(K)  12.6465 

CONDUCITVITY  ,3880 

PRANDTL  NUMBER  12.1432 

WITH  FROZEN  REACTIONS 

CP,  CAL/(G)(K)  .3514  <— 

CONDUCITVITY  .1781 

PRANDTL  NUMBER  .7349 


The  mole  fractions  for  this  calculation  are  the  same  as  those  in  Appendix  A  and  are  not  reiMDduced  here. 
The  amount  of  graphite  produced  is  small  enough  to  warrant  ruiming  this  calculation  with  the  NASA- 
Lewis  code.  Cy  is  calculated  according  to  Equation  18  using  the  Cp  returned  in  the  transport  properties 
calculation  with  frozen  reactimis  (.3514  cal/g-K),  M  (29.013  ^hrol),  and  R  (1.987  cal^sol-K).  AT  is  the 
difference  between  T(^  and  298.15  K: 

Cy  -  .3514  -  1.987/29.013  -  0.2829  cal/g-K 
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AT  -  900  -  289.15  -  601.85  K 
CvAT  -  170.27  cal/g 

The  calculated  internal  energy  of  the  combustion  products  is  the  difference  between  U  calculated  at  Tfo 
and  CyAT: 

^298. 15  “  “  1259.77  -  170.27  -  -  1430.04  cal/g 

The  internal  energy  of  the  imcombusted  HMX  is  subtracted  from  this  value  to  determine  HEX.  A 
correction  of  40  cal/g  is  added  to  this  value  for  the  CHNO  chemistries  looked  at  in  this  report: 

HEX  =  1-1430.04  -  +60.511  =  1491  cal/g 
HEX  +  40  =  1491  +  40  =  1531  cal/g 
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APPENDIX  C: 

EXAMPLE  OF  ROOM  TEMPERATURE  METHOD  FOR  HMX  AT  A  LOADING  DENSITY 

OF  0.1  g/cm^ 
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Prior  to  running  this  method,  the  thenno  library  could  be  restricted  to  CO.  CO2,  H2.  N2.  H2O,  and 
H20(1)  to  eliminate  the  prediction  of  unrealistically  large  amounts  of  C(gr),  CH4,  and  NH3.  Run  an 
assigned  temperature  and  volume  (TV  =  T)  calculation  at  room  temperature.  The  input  file  looks  as 
follows: 


REACTANTS 

C  4.  H  8.  N  8.  08.  .  1.0  17924. 5S  298.15  HMX 

NAMELISTS 

&INPT2  TV-T,T-298.15,RHO-0.1,TRACE=1.E-03/ 


The  following  calculated  values  are  extracted  from  the  output  file: 


THERMODYNAMIC  PROPERTIES 


P,  ATM 

T,  DEG  K 
RHO,  G/CC 
H,  CAL/G 

U,  CAL/G 
G,  CAL/G 

S,  CAL/ (G) (K) 

M,  MOL  WT 
{DLV/DLP)T 
(DLV/DLT)P 
CP,  CAL/(G)(K) 
GAMMA  (S) 

SON  VEL, M/SEC 
IMPETUS  (J/G) 


MOLE  FRACTIONS 

CO 

C02 

H2 

H20 

N2 

H20(L) 


98.833 
298.15 
1.0000-1 
-1270.23 
-1294.17  <— 

-1702.13 
1.4486 

24.755 

-1.00361 

1.0094 

.3166 

1.3419 

366.6 

100.1 


3.3352-3 

3.3000-1 

3.3000-1 

.00031 

3.3333-1 

3.0218-3 


Since  the  calculation  is  tun  at  room  temperature,  it  is  a  simple  matter  to  subtract  the  internal  energy  of 
the  uncombusted  HMX  from  the  calculated  U  to  get  HEX: 


HEX  -  1-1294.17  -  +60.511  -  1355  cal/g 
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Intentionally  left  blank. 


30 


APPENDIX  D: 

EXAMPLE  OF  THE  INERT  DILUENT  METHOD  FOR  HMX 
AT  A  LOADING  DENSITY  OF  0.1  g/cm’ 


31 


IhrTENTIONALLY  LEFT  BLANK. 
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Prior  to  running  this  method,  the  thermo  library  could  be  restricted  to  CO.  CO2.  H2.  N2.  H2O.  H20(l). 
and  Argon.  Run  an  assigned  energy  and  volume  (UV  =  T)  calculation  ^cifying  the  energetic  material 
(or  formulation)  and  Argon.  The  following  example  will  use  1  part  HMX  to  999  parts  Argon  by  weight 
The  input  file  looks  as  follows: 


REACTANTS 

C  4.  H  8.  N  8.  08.  .  1.0  17924. 5S  298.15  HMX 

ARl.  ....  999.0  -592. 5G  298.15 

NAMELISTS 

&INPT2  UV-T,RHO-0.1, TRACE-1. E-07/ 

The  following  calculated  values  are  extracted  from  the  output  file: 


THERMODYNAMIC  PROPERTIES 


P,  ATM 

T,  DEG  K 
RHO,  G/CC 
H,  CAL/G 

U,  CAL/G 
G,  CAL/G 

S,  CAL/(G)(K) 

M,  MOL  WT 
(DLV/DLP)T 
(DLV/DLT) P 
CP,  CAL/(G)(K) 
GAMMA  (S) 

SON  VEL, M/SEC 
IMPETUS  (J/G) 


MOLE  FRACTIONS 

AR 

CO 

C02 

H2 

H20 

N2 


65.003 
316.26  <-= 

1.0000-1 
.98551 
-14.756 
-228.88 
.7268 

39.923 
-1.00000 
1.0000 
.1246  <— 

1.6656  <— 

331.2 
65.9 


9.9838-1 

2.6729-6 

5.3655-4 

5.3655-4 

2.6729-6 

5.3922-4 


HEX  is  calculated  directly  from  Equation  20  (N  ==  KXX)): 

HEX  -  1000*0. 1246*(316. 26-298.15)/!. 6656  -  1355  cal/g 
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APPENDIX  E: 

DERIVATION  OF  EQUATIONS  FOR  BOMB  CALORIMETER  COMBUSTION 
AND  EXPLOSION  REACTIONS 
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iNTEhrnONALLY  UEFT  BLANK. 
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This  Appendix  details  how  the  molecular  coefficients  in  Equations  8  and  10  in  the  body  of  this  report 
were  (Stained.  Combustion  of  a  compound  in  an  excess  of  O2  would  be  expected  to  go  to 

completion  and  yield  CO2,  H2O,  and  N2,  along  with  any  umeacted  02-  For  this  combustion  reaction,  and 
neglecting  the  excess  O2.  start  with  the  following  ui^anced  chemical  equation: 


C.HbOcNd  +  WO2  ->  xCOj  +  yH20  +  ZN2.  (E-1) 

The  coefficients  x,  y.  and  z  in  Equation  E-1  are  easily  determined  from  tire  a,  b.  and  d  coefficients 
since  C.  H,  and  N  each  appear  once  on  the  right  side  of  the  equation:  z  =  d/2,  y  =  b/2.  and  x  =  a. 
Equation  E-1  now  becomes: 

C.HbOcNj  +  WO2  aC02  +  IH2O  +  £n2.  (E-2) 

The  O2  coefficient,  w,  can  now  be  solved  for  in  terms  of  the  coeflicients  for  those  species  which  contain 
oxygen.  If  the  coefficents  are  put  in  terms  of  O2.  or  two  atomic  oxygen  equivalents,  the  following 
equation  can  be  used  to  solve  for  w: 

£  ♦  w  -  a  ♦  (E-3) 

2  4 

where  c/2  comes  from  C,H|,OcNj,  w  is  the  coefficient  for  O2.  a  is  the  coefficient  for  CO2,  and  b/4  is  one- 
half  the  coefficient  for  H2O.  Equation  E-3  is  easily  solved  for  w.  Making  the  substitution  for  w  in 
Equation  E-2  yields  Equation  8: 

C,Hb0^d  +  (a  +  ^  -  ^)02  aC02  *  JH2O  *  |n2.  (E-4) 


Equation  10  can  be  derived  in  a  similar  manner.  The  explosion  process  here  is  treated  as  an 
incomplete  combustion,  since  the  reaction  proceeds  under  an  inert  environment  such  as  Ar  or  N2.  In  this 
case,  only  H2O,  CO2.  N2.  and  CO,  will  be  considered  as  combustion  products.  The  explosion  reaction 
can  now  be  written  in  an  unbalanced  form  as: 

C,Hj,OcN,j  -»  WCO2  ♦  xCO  ♦  yH20  +  ZN2.  (E-5) 
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The  H2O  and  N2  coefficients  are  determined  in  exactly  the  same  way  as  was  the  case  for  the  combustion 
reaction  above:  y  =  h/2.  z  =  d/2.  This  leaves  the  CO2  and  CO  coefficients,  w  and  x,  which  can  be 
determined  from  the  a,  b,  and  c  coefficients.  The  partially  balanced  reaction  looks  like 

C.HbO^Nj  -»  WCO2  +  +  3H2O  +  1n2.  (E-6) 

The  following  two  equations  are  set  up  to  conserve  the  C  and  O  mass  balance; 

C:  a  -  w  +  X  (E-7) 


O: 


c  =  2w  +  X  +  _, 
2 


(E-8) 


where  the  coefficients  a  and  c  in  Equation  E-7  and  E-8  come  from  the  w,  x.  and  b  coefficients 

come  from  CO2.  CO,  and  H2O,  respectively,  and  the  O  mass  balance  in  Equation  E-8  is  based  on  one 
atomic  oxygen  equivalent  One  can  solve  for  w  by  subtracting  Equation  E-7  from  Equation  E-8.  It  is 
then  a  simple  matter  to  make  the  substitution  for  w  in  Equation  E-7  to  solve  for  x: 


w  *  c  -  a 


b 

2 


(E-9) 


X 


-c  +  2a  ♦  -. 

2 


(E-10) 


Making  the  substitutions  for  w  and  x  in  Equation  E-6  yields  Equation  10: 

C.HbO^Nd  -»  (c  -  a  -  ^)C02  +  (-c  +  2a  +  .^)CO  +  J^H20  +  £n2.  (E-H) 

2  2  2  2 
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